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ABSTRACT 

We show how the redshift and peak-flux distributions of gamma-ray bursts (GRBs) 
have an observation-time dependence that can be used to discriminate between dif- 
ferent burst populations. We demonstrate how observation-time relations can be de- 
rived from the standard integral distributions and that they can differentiate between 
GRB populations detected by both the BATSE and Swift satellites. Using Swift data 
we show that a redshift-observation-time relation (log Z - log T) is consistent with 
both a peak-flux - observation-time relation (log P- log T) and a standard logiV- 
logP brightness distribution. As the method depends only on rarer small- z events, 
it is invariant to high-z selection effects. We use the log Z- log T relation to show 
that sub-luminous GRBs are a distinct population occurring at a higher rate of or- 
der 150^9® Gpc yr -1 . Our analysis suggests that GRB 060505 - a relatively nearby 
GRB observed without any associated supernova - is consistent with a sub-luminous 
population of bursts. Finally, we show that our relations can be used as a consistency 
test for some of the proposed GRB spectral energy correlations. 

Key words: gamma-rays: bursts - gamma-ray: observations - methods:data analysis 
- supernovae: general - cosmology: miscellaneous 



1 INTRODUCTION 

Multi-wavelength observations have shown that 7-ray bursts 
(GRBs) are the most luminous ]] and distant transient event s 
, in the Universe (|Greiner et al. 1200 i lCucchiara et al.ll201lh . 
GRBs have been generally categorized into two populations: 
spectrally soft long duration bursts related to core-collapse 
events (LGRBs/Type II); harder short duration bursts pos- 
sibly resulting from compact star mergers (SGRB/Type I). 

In addition to these two main populations of bursts 
it has been suggested there exist two sub-populations: 
sub-luminous GRBs (SL-GRBs) and SGRBs with ex- 
tended emissions (SGRB-EE). SL-GRBs are of the long- 
duration type and have isotropic equivalent 7-ray lumi- 
nosities 2-3 or d ers of magnitude below classical LGRBs 
JCowardl [20051: iMurase et all 120061 : iGuetta fc Delia Valid 
l2007l : llmerito et al.ll20Qg| ). The lower energy emissions mean 
they are only detected at low- z - as such, four of the six 
GRBs with unambiguous spectroscopically confirmed GRB- 
supernova associations were from this category. SGRB-EEs 
emissions have been gi ven a separate cla ssification in the sec- 
ond Swift catalogue (|Sakamoto et al.l 120111 ). These bursts 
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In terms of electromagnetic radiation per unit solid angle. 



show an initial SGRB like short hard s pike ( < 2s) fol- 
lowed by a faint softer emission (> ,100 s) dNorris fc Bonnelll 
| 200d:lNorris et al.ll201ll : iPage et alj|200a iPerlev et al.ll2008l : 
Izhandl20lj )~ 

There is still no clear consensus that these sub- 
categories arise from different progenitor systems or are sim- 
ply rarer events from the tail of the respective short/long 
burst distributions. Attempts to add ress this have generally 
been based on statistical ar guments l|Soderberg et al]|2006l : 
IGuetta fc Delia Valid 120071 ) . fits to the l og N- log P. peak 
flux, or 'brightness dist ribution ' of bursts l|Pian et al] 20061) 
or through simulation l|Cowardll2005l : IVirgilii et al.ll20od h 

The goal of this paper is to demonstrate an alternative 
strategy using the relative time records of the bursts. This 
approach exploits the fact that different astrophysical tran- 
sient populations will have different local rate densities. We 
show that by recording the arrival times of the rarest events 
in a time series e.g. the closest or brightest of a cosmological 
population, one can produce a rate dependent data set with 
a unique statistical signature (jCoward fc Burmanll2005h . By 
constraining the data using an observation-time dependent 
model that is highly sensitive to the rate density, we demon- 
strate how this alternative approach can untangle different 
source populations. In this study we will specifically address 
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the issue of whether SL-GRBs are a distinct population of 
GRBs. The outline of the paper is as follows: 

In section 2 we present an overview of GRB popula- 
tion studies. Section 3 will set the scene in regards to the 
observation-time dependence of transient events and Sec- 
tion 4 will describe the data extraction methodologies used 
in this study. A standard theoretical framework will follow 
in Section 5. 

Section 6 describes observation-time dependent models 
for both peak flux (log P- log T) and redshift (log Z- log T) 
showing how they follow seamlessly from the relative integral 
distributions of transient sources. In Section 7, parameters 
for both the BATSE0 and Swift LGRB populations will be 
obtained using a standard differential log N— log P distribu- 
tion. These parameters will be used in Section 8 to constrain 
the peak-flux - observation-time data from both detectors 
using a log N- log P model. Doing so demonstrates that the 
method is both consistent with a standard brightness distri- 
bution and is detector independent. 

Section 9 will demonstrate the use of the previously 
derived parameters in the redshift - observation-time do- 
main to constrain Swift data using the log Z- log T model. 
In Section 10 we apply our methods to the Swift SL-GRB 
population to further demonstrate how the method can be 
used as a tool to discriminate between different source pop- 
ulations. We show that the method uses only the closest or 
brightest of a population; thus many of the selection biases 
that plague GRB observations can be bypassed. We discuss 
our findings and present our conclusions in Section 11. 



2 GAMMA-RAY BURST POPULATIONS 

The catagorisation of GRBs was traditionally based on the 
bi modal distribution of Tg p durations observed by BATSE 
(|Kouveliotou et all 1 1993! ) and their hardness ratio in the 
spectral domain. These criteria separated GRBs into hard 
SGRBs (T 90 < 2 s; hard spectra) and softer LGRBs (T 90 ^ 
2 s; soft spectra). 

Electromagnetic observations of LGRBs and SGRBs 
have provided strong evidence for different progenitors. 
LGR Bs have been associated with the dea t hs of massive 
stars l|Wooslev fc Blooml |2006l : iHiorthl 120031 ; IStanekl [2003 ) 
and have subsequently been found in or near dense regions 
of active star-formation, predo minantly dwarf starburst field 
galaxies |Fruchter et alj |2006). For SGRBs, which have con- 
tributed around 25% and 10% of the BATSE and Swift 
GRB samples respectively (|Guetta fc St ella 2009), the lead- 
ing progenitor model is the merger of compact neutron stars 
and/or black hole binaries. The association of an older stel- 
lar population with these bursts is supported by their oc- 
currence in both early and late-type galaxies, as well as field 
and cluster galaxies. 

There exist, however, a number of ambiguities in the 
categorisation schemes of GRB populations. LGRBs, such 



2 The Burst and Transient Source Experiment (BATSE) on the 
Compton Gamma-Ray Observatory was launched in 1991, and 
recorded over 2000 GRBs at a rate of around one a day for around 
9 years of operation. 

3 The Tgo duration is the time during which the cumulative 
counts increase from 5% to 95% above background. 



as GRB 060614 and GRB 060505, showed no evidence 
of a supern ova, despite extensive follow up campaigns 
(|Zhangir2006l '). Additionally, it has been suggested that the 
two closest recorded bursts, GRB 980425 (36 Mpc) and 
GRB 060218 (145 Mpc), along with GRB 31203, asso- 
ciated with a host galax y at ~ 480 Mpc ( Feng fc Foxl 



l2010h and GRB 100316D l|Starling et all l201ll)fz ~ 0.06 



make up a sub-class of SL-GRBs ( Cobb et al.l I200 1 



Liang et all 120071 : iGuetta fc Delia Valid 



20071 ; IVirgili et all 



2009). This class of GRB have isotropic equivalent 7-ray en- 



ergy emissions typically several orders of magnitude below 
those of standard long-duration GRBs (jMurase et al.l 120061 ; 
IGuetta fc Delia Valldl2007l ; llmerito et al l l2008h suggesting 
that they could form a unique population of bursts that due 
to their relatively close proximity, must be occurring at a 
higher rate. The suggestion that SL-GRBs could be just nor- 
mal LGRBs viewed off-axes has also been considere d . This 
scenario was discounted by iDaigne fc Mochkovitchl (|2007h 
on statistical arguments as it would: a) result in a far higher 
local rate density than expected from LGRBs ; b) require 
much narrower opening angles for LGRBs than typically de- 
rived from the breaks in afterglow lightcurves. 

While the LGRB/supernovae connection is firm for the 
majority of these bursts, the progenitors behind SGRBs, 
which are rarer and more difficult to localise, are less cer- 
tain (|Meszaroj 120061 : iNakaij l2007l i. Around 20% of the 
SGRBs detected by Swift have b een followed by an ex- 
tended emission lasting up to 100 s jNorris fc Bonnelll f2006: 
IPerlev et al.ll200ljl ) leading to su ggestions that diffe rent pro- 
genitors produce these bursts (|Norris et alj|201lh . Candi- 
date systems include the birth of a rapidly rotating proto- 
magnetar produced through NS-NS merger or the accre- 
tion induced collapse o f a white-dwarf |Metzger et al.|[2008l : 
iBucciantini et aTll2012h . IChapman et al.l (|2009l ) further pos- 
tulated that an initial spike produced by a soft 7-ray re- 
peater (SGR) in a galaxy of close proximity could mimic a 
SGRB and found that a dual population luminosity func- 
tion based on both SGR giant flare properties and SGRB 
luminosities was co nsistent with BATSE data. Interestingly, 
IVavrek et al.l ((2008) found the sky distribution of SGRBs to 
be anisotropic, implying a fraction at close proximity. Cer- 
tainly, this suggests that classificati on of shorter du ration 
bursts could be detector dependent (|Page et al.ll2006l l. 

The various ambiguities have motivated a number of 
authors to re-define different classes of GRB through a 
number of properties including: spectral features, associated 
supernova, stellar population , host galaxy, location in the 
host galaxy and progenitor (|Zhang et al ] |2007I :IbI 00m et al.l 
2008). This has led to the Type I (compact object mergers) 
and Type II (core-collapse) scheme. 

The existence of a third population of bursts with inter- 
mediate duration has also been suggested. The first evidence 
of this additional population was provid ed though statis- 
tical analysis of the BATSE distribution (M ukheriee et all 
ll99StlHorv ath 1998). This was la ter supported b y analysis 
using data f rom t he BeppoSAX (|Horvkth|[20Qgh. RHESSI 
(Ri'pa ~et al.l |2009h and Swift satellites ( Horvath et al.l 
l2010al ). To for mulate a phy s ical m odel for an intermediate 
class of bursts, IVeres et al.l (|201CJ ) showed that intermedi- 
ate bursts have a lower than average peak flux distribution 
and suggested that this group, which are spectrally softer 
than LGRBs, may be related to X-ray flashs (XRFs) - these 
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events have a dominant fraction of the total prompt fiuence 
detected in X-ray rather than 7-ray. 



3 EXPLOITING THE TIME DIMENSION 

In this study we show that the time-record of GRB observa- 
tions can be used as a tool to untangle different GRB pop- 
ulations. This largely overlooked quantity features strongly 
when one considers the observation of the rarest events in 
a population i.e. those events from the tail of the distri- 
bution which occur at low-2 or have exceptional brightnes s 
in comparison with the average. ICoward fc Burmani (|2005l ) 
showed that these events posses a unique rate dependent 
statistical signature that can be described by the 'probabil- 
ity event horizon' (PEH) concept. By recording successively 
rarer events as a function of observation-time, a data set 
can be produced and constrained by a rat e dependent model 
i|Coward et al.ll2005l ; iHowell et alll2007al ). 

The basic concept is illustrated in Figure [1] showing the 
Swift LGRB peak flux data plotted against observation-time 
- we will often refer to such a time-series as P(T) data (or 
Z(T) in the redshift domain). Using a log- log plot it is ap- 
parent that successively brighter events have an observation- 
time dependence - the longer you observe, the greater the 
probability of observing an exceptionally bright event. Suc- 
cessively brighter events are indicated by squares - these are 
termed PEH data. 

To constrain PEH data lHowell et all (|2007bl) extended 
the standard Euclidean derivation of a log N- log P bright- 
ness distribution to show how including a time dependence 
could produce a peak flux-observation-time power law rela- 
tion, log P- log T. They also illustrated how some processing 
could improve the amount of PEH data. In an application to 
Swift LGRB data, they demonstrated that the relation could 
be used as a probe of the event rate density of the sources. 
Additionally, the time dependence allowed the method to be 
used as a tool to predict the likelihood of future events. 

In this paper we extend the previous study to present a 




observation time (days) 

Figure 1. The Swift LGRB peak flux data plotted against 
observation-time. The plot illustrates how the probability of ob- 
serving a bright event increases with observation-time. Succes- 
sively brighter events - termed PEH data - are indicated by 
squares. 



cosmological log P- log T relation and show that the princi- 
ple also extends to the redshift domain resulting in a log Z- 
log T relation. For the latter case, the PEH data is a measure 
of the geometrical distribution of the source population. A 
log Z— log T curve can be thought of as a horizon, defined 
by successively closer events, that app roach a central ob- 
serve r as a function of observation-time (jCoward &l BurmarJ 
l2005h . As the successively closer events approach the local 
low- 2 : regime rapidly, the GRB selection function l|Cowardl 
2007) and high-z sele ction effects such as the 'redshift desert' 
( Coward et al.ll200ct ) have a negligible effect. 

In the next section we will describe the data extraction 
principles that will be employed in the latter part of this 
study. 



4 DATA EXTRACTION METHODOLOGY 

To extract PEH data, w e follow the FromMax method used 
bv lHowell et all (|2007bl ) to constrain the Swift LGRB sam- 
ple. This invoked the temporal cosmological principle: for 
time scales that are short compared to the age of the Uni- 
verse, there is nothing special about the time we switch on 
our detector. Therefore, the P(T) time-series can be treated 
as closed loop, i.e. the start and endpoints of the time se- 
ries can be joined and the start time set immediately after 
the brightest event. Successively brighter events and their 
observation-times are then recorded to produce a PEH data 
set of peak fluxes. 

Figure EQ illustra tes the concept. Using this procedure 
IHowell et afl l |2007bl ) showed that PEH data set can be ob- 
tained through both temporal directions - for simplicity we 
will only record data in the forward direction. 

Employing this technique circumvents the possibility of 
a loud event occurring early in a time series; this would min- 
imize the amount of output data as the next largest event 
would most probably occur near the end of the time series. 
Such a situation could be encountered by a detector operat- 
ing with a high energy cutoff - a bias could be introduced 
producing a large number of events reaching the threshold. 
The advantage will become apparent in section I5TT1 through 
analysis of the BATSE sample of bursts which contained a 
number of bright events at early observation-times. 

Another feature of the FromMax method is that it en- 
sures the total time duration of the PEH output is always 
equivalent to that of the total observation-time - this en- 
sures a well ordered data sample is produced with a con- 
sistent time signature. We will show later through 2D Kol- 
mogorov Smirnov testing that the improved data set retains 
the statistical signature of the original. 

To apply the procedure to a sample of P(T) data we 
first define the brightest event by P* with an observation- 
time stamp T* and denote the time of the last, most recent 
occurring event, as Tmax- Treating the data as a closed loop 
we reorder the data starting from the first event after P*. 
The time stamps of the re-ordered data set P\Tp) are now 
defined as: 

, T - T» T >T, 

T'p = (1) 
L T + T max — T» T^T, 



To obtain a PEH data set we start to extract our data 
from the first minimum P^in = Pi < Pl+i ~ this additional 
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Figure 2. A cartoon illus trating the FromMAX data extraction 
principles of lHowell et al,l 1 1200711 ). Treating the data as a closed 
loop one can select a start point from which a data set of suc- 
cessively brighter events can be extracted. This data set is called 
PEH data. 



step is to minimise the effect of an early bright event. A PEH 
data set is then obtained by recording successively brighter 
events (P|, Tp i ) satisfying the condition Pj+i > P for Pi ^ 
P" ■ 

' min ■ 

To determine the PEH data set in the redshift domain 
one applies similar principles, treating the data as a closed 
loop but re-ordering the data from the first event after the 
closest redshift event Zq. The timestamps for the re-ordered 
set Z K (T Z ) are then given by: 



T-To 



T > To 



(2) 



T + T max — To T ^ To 



A PEH data set is obtained by extracting data from the 
first maximum Z x max = Z\ > Z\ +1 , recording successively 
closer events (Z},Tz A ) satisfying the condition Z\ +1 < Zi K 
for Z; 5? Z m ~^. 



5 THEORETICAL FRAMEWORK 

5.1 GRB Flux and luminosity relations 

We firstly define an isotropic equivalent photon luminosity 
in the source frame as: 

/■10000 kcV 

L= S(E)dE, (3) 

AkcV 

where S(E) is the differential rest-frame source luminos- 
ity (in ph s _1 keV -1 ). To define the observed peak pho- 
ton flux (photons per cm 2 per second) observed within 
a detector band E m i n < E < -E max and emitted by an 
isotropically radiating source at redshift z one must 
perform two modifications. Firstly, the observed photon 
flux is modified to account for the missing fraction of 
the gamma ray energy seen in the detector band b = 



Iei S(E)dE/ J* 0000 S(E)dE. Secondly a cosmological k- 
correction, k(z) = S{E)dE/ f^*"} S(E)dE is ap- 
plied. With these two modifications the standard definition 
for flux becomes: 

L b 



(l + *0 



47tc(l(z) 2 k(z) 



(4) 



Substituting in for b and k(z) one obtains the familiar rela- 
tion: 



P = 



(! + *)/(' 



(l + 2)B max 

l+z)B min 



S{E)dE 



A-Kd^z) 2 



(•>) 



where di J (z) is the luminosity distance. The (1 + z) fac- 
tor is included as the standard definition of dh{z) is valid 
for an energy flux, b ut P here is given as a photon flux 
l|Meszaros et al.ll201ll ). For long duration GRBs the func- 
tion S(E) is typically modeled by a Band function l|Bandl 
which we use with high and low energy spectral in- 
dices of -2.25 and -1 and a break energy of 511 keV. 



5.2 GRB Luminosity Function 

A number of different forms have been suggested for 
the LGRB Luminosity Function (LF). To minimise our 
free parameters we the single power law form used by 
IPorciani fc Madaul l|2000l ) which has an exponential cutoff 
at low luminosity: 



<]>(/.) = <!>„ ( ^~ 



exp - — 



(6) 



Here, L is the peak rest frame photon luminosity in 
the 1-10000 keV energy range, a ensures an asymptotic 
slope at the bright end and L, is a characteristic cut- 
off scaling. The normalisation coefficient for this LF is 
gi ven by $p = \L*T(g — 1 ) ]~ . B ased on the studies 
oflMeszaros fc Meszarosl dl99dll995l '); lRe'ichart fc Meszaroj 
\ 19971 ): iButler et al. I (|201oF we will assume no luminosity 
evolution with redshift. 



5.3 GRB rate evolution 

To model GRB rate evolution, -Rgrb(z), we assume that 
LGRBs track the sta r formation history of the Universe 
|Meszaros et al.ll2006l 1 and normalise a star formation rate 
model, Rsf{z), (in units of mass converted to stars per unit 
time and volume) to the local (z = 0) rate. Multiplying 
by the local rate density po (in Gpc _3 yr _1 ) allows one to 
extrapolate rate evolution to cosmological volumes: 

Rsf(z) 



Rgrb(z) = po- 



(7) 



Rsf(z = 0) ' 

For Rsf(z) we use the model of|lj| l|2008f l. which was ob- 
tained by addi ng ultraviolet and i n fared measurements to 
the sample of iHopkins fc Beacoml (|2006t ). This takes the 
form: 

(0.02 + 0.122) 



Rsf(z) = 



1 + (z/3.23) 4 - 66 ) 



A number of studies suggested that cosmic metallic- 
ity evolution must also be considered in any rate ev olu- 
tion model for LGRBs (jll l200g| ; iModiaz et all 120081 ). A 
metallicity dependence results from the requirement that 
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Wolf-Rayet (WR) stars should retain sufficient rotation 
to power a GRB - therefore angular momentum losses 
through stellar-wind ind uced mass-loss must be minimized 
l|Wooslev 8z Heg cr 2006). As wind-driven mass loss in WR 
stars is understood to be dependent on a high enough frac- 
tion of iron , a low-metallicity environment is an essential re - 
quirement (jvink fc de Koterll2005l; IWooslev fc Jankall2005l) . 

Recent results from lElliott et al. l l|2012i ). driven by ob- 
servations of LGRBs in metal rich galaxies, suggest that in- 
terpretations of luminosity and metallicity evolution could 
simply result from a misu nderstanding of various redshift 
biases l|Coward et al]|2~008l ). As we will demonstrate later in 
section [9l the methods employed in this paper are sensitive 
to the rarer low- z population of bursts. Therefore, for clar- 
ity, and in support of the uncertainties discussed above, we 
will ignore the effects of metallicity evolution in this study. 




1 2 3 4 5 

Log T (days) 

Figure 3. Cosmological and Euclidean log P- log T curves. The 
two curves converge at high observation-times. 



5.4 The all sky event rate equation of GRBs 

The number of GRBs per unit time within the redshift shell 
z to z + dz with luminosity L to L + dL is given by: 

dA 



dtdzdL 



Here the (1 + z) factor accounts for the time dilation of the 
observed rate by cosmic expansion; its inclusion converts 
source-count information to an event rate. The co-moving 
volume element: 

dV _ 4nc d L 2 (z) 
dz 



(10) 



Ho (l + z)2h(z) ' 

describes how the number densities of non-evolving objects 
locked into Hubble flow are constant with redshift. The 
quantity h(z), is the normalized Hubble parameter, 



h(z) = H{z)/H = [fi m (l + zf + fi/ 



,1/2 



(11) 



wher e f2 m + Qa = 1 (for further details see ICarroll et al.l 
1 1992! ) . For a 'flat-A' cosmology, we take fi m = 0.3, f2 A = 0.7 
and Hq = 72 km s _1 Mpc -1 for the Hubble parameter at 
the present epoch. 



6 THE OBSERVATION-TIME RELATION FOR 
PEAK FLUX AND REDSHIFT 

6.1 The log P— log T relation 

From equation [5] the rate of GRBs with a peak photon flux 
greater than P observed by an instrument with sky coverage 
is given by: 



AL,P) 



dV(z) R G rb(z) 
dz (1 + z) 



dz, (12) 



were z ma x(£, P) is the maximum redshift from which a burst 
with luminosity L and peak flux P can be detected. 

To introduce an observation-time dependence, T, 
we follow the p r obabi lity event horizon concept of 
ICoward fc Burmanl (|2005T ) and note that as GRBs are in- 
dependent of each other their observation-times will follow 
a Poisson distribution in time. Therefore, the temporal sepa- 
ration between events will follow an exponential distribution 



defined by a mean number of events, AT(> P) T. The prob- 
ability e for at least one event > P is given by: 



P(n>l;N(>P),T) 



N(>P)T 



(13) 



For this equation to remain satisfied with increasing 
obser vat ion- 1 ime : 



N(> P)T= |ln(l — e)| 



(14) 



Equating the above equation for P and T we obtain a re- 
lation for the evolution of peak flux as a function of obser- 
vationtime. By setting e to some arbitrary value, log P- log 
T curves can be obtained numerically through equations 1121 
andlHl 

Assumin g a con stant radial distribution of sources, 
iHowell et all l|2007bl ) derived the following compact form 
for a log P- log T relation, which is a T 2 / 3 power law for a 
Euclidean distribution of sources: 



P(T)= T 



2/3 



Aflro 



3V47r|ln(l - e) 



2/3 



${L)L i/2 dL 



2/3 



(15) 



Figure [3] compares both the cosmological (equations [12] 
and !14[) and Euclidean fequation ll5|) curves for arbitrary pa- 
ramete rs of source evolution and LF. Following lHowell et al.l 
( 2007a|) we have plotted upper and lower thresholds by set- 
ting e = (0.95; 0.05) - we will refer to these curves as 90% 
PEH thresholds. We note that as observation-time increases 
the Euclidean and Cosmological curves begin to converge. 
This represents the increased probability of obtaining a cos- 
mological bright event from the bright —3/2 power law part 
of the log A-logP distribution. The Euclidean curves are 
convenient in that they can provide initial estimates by fit- 
ting to the brightest sample of PEH peak flux data. 



6.2 The log Z— log T relation 

One can extend the arguments of the previous section to 
derive a log Z- log T relation. From equation [5] the rate of 
GRBs observed by an instrument with sky coverage Q. within 
a redshift limit is given by: 
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imi„(PL,Z L ) 



dV(z) Rgrb(z) 
dz (1 + z) 



dz , 



(16) 



with obtained by applying the flux limit of the detector, 
Pl, to equation [S] The quantity 7? z is the efficiency of ob- 
taining a redshift, approximated as the fraction of the total 
burst sample with measured redshifts. Using a similar ar- 
gument as used to determine equation 1141 one obtains the 
following relation for the temporal evolution of redshift: 



N(< Z L )T = |In(l-e)|. 



(17) 



This equation can be equated for T and z to set a spatial 
dependence on GRB populations. Curves of log Z— log T for 
e = (0.95; 0.05) can be obtained numerically through equa- 
tions \W\ and 1171 

By using a single lo g Z-logT curve as a threshold 
ICoward fe Burmanl ( 20051 ) showed that SL-GRBs appeared 
as outliers, supporting the suggestion that they were a sub- 
population of classical LGRBs. In this paper we have ex- 
tended their relation, based solely on the cosmological event 
rate evolution Rgrb(z), to include the effects of detector 
sensitivity and the luminosity distribution of sources. As 
shown above, this allows the log Z-logT relation to be de- 
rived seamlessly from a standard integral distribution. Thus, 
model parameters obtained by fitting to a differential bright- 
ness distribution should satisfy the two observation-time 
relations presented above (equations [14] and I17|l . This cir- 
cumvents a difficulty encountered in previous studies which 
fitted directly to the PEH data to determine a rate - the 
resulting rate estimates h ad low resolution due to the intrin- 
sic sc atter of PEH data (|Howell et al.l l2007al : iHowell et ail 
l2010l) . 

Before we go on to apply the fits described in sec- 
tions 16.11 and 16.21 to PEH data, we will first determine 
the parameter values that will be used in the models of 
log (P; Z)-log T. We will do this by fitting our data to the 
brightness distribution of GRBs. 



7 FITTING TO THE GRB BRIGHTNESS 
DISTRIBUTION 

In this section we perform \ 2 minimisation of the bright- 
ness distribution of both Swift and BATSE LGRBs to ob- 
tain values for the three free parameters po, Lo and a. A 
goodness of fit for this procedure is obtained by dividing 
the minimised \ 2 value by the number of degrees of free- 
dom, xV dof - We note that although the log (P; Z)- log T 
relations are derived from an integral distribution, we 
fit to a differential distribution in which the number of 
sources are independent at each interval of P + dPQ. 
From equ ation [9] one can define a differential log A'- log P 
relation feommers et all 120001; iPorciani fc Madaul |2000| ; 
ISalvaterra fc Chincarinill2007l ; ICampisi et al.ll2010l ) which is 



4 Bright objects will contribute to counts at all values of P in an 
integral distribution 




log P ( ph cm " s 1 ) 

P distribution of BATSE 



Figure 4. The differential log N 
(top) and Swift (bottom). 



the observed rate of bursts within a peak flux interval (Pi , 
Pi) as: 

£(-Pl,z) 

with L(Pi i 2,«) obtained though equations 2] and \S\ 

We bin peak flux data into logarithmical spaced in- 
tervals AP and ensure each bin contains at least 5 bursts 
l|Wall et alj|2003ft . Bursts per bin AN and their uncertain- 
ties ±V A7V are converted into bu rst rates AP by divid ing 
by the live time of the search AT (|Kommers et al.M2000f ). 



7.1 The BATSE log TV - log P distribution 



A nu mber of studies have utilized the iKommers et al.l 
2000) differentia l peak flux distribution dat a (eg . 



Porciani fc Madau 2000; ISalvaterra fc Chincarinil [2007; 



Campisi et al.l 2010). As we require both the discrete peak 



flux values and their observation-times to produce a P(T) 
distribution in the next section, we make our own data se- 
lection fr om the BATSE c urrent catalogu^S We note how- 
ever, the IKommers et al.l (|2000h sample is a good consis- 
tency test for our models and we find good agreement 
using the best fit parameters of (|Porciani fc Madaul l2000l : 
ISalvaterra fc Chincarinill2007l ). 

The BATSE instrument performed a number of runs 
at different trigger energy ranges and photon count rate 
thresholds throughout its operation. These runs, in effect, 
represented different experiment^. For consistency, we con- 
dition this data by selecting bursts from the (50-300) keV 
band obta ined when the trigge r thr esholds were set to 5.5g . 
Following iGuetta et all l|2005h and iGuetta fc Piranl (|2007n 
we additionally select bursts for which the ratio of count 
rates, C, at the 1024 ms timescale Cmax/Cmin ^ 1 — here 



5 http : //www. bat se .msfc .nasa.gov/batse/grb/catalog/current/ 

6 For example, during the first 17 months of operation the energy 
range was set to (50-300)keV, although the trigger thresholds were 
varied. 
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Data 


Trigger 
numbers 


Number of 
bursts 


Run time 
(days) 


A 


105-178 


11 


18 


13 


268-1851 


164 


437 


C 


1928-3175 


198 


724 


D 


3883-3941 


14 


41 


E 


5403-5519 


37 


76 


F 


6102-6764 


143 


448 


G 


7356-7767 


29 


228 



Table 1. The BATSE data segments of bursts that triggered in 
the (50—300) keV band when the trigger thresholds were set to 
5.5cr. We show the number of bursts within each segment that 
passed our selection criteria and the instrument run time associ- 
ated with each segment. 



20071) . IWanderman fc Piranl (|201Ch and Salvaterra et al.l 
2009) and adopt a simplified approach to account for de- 
tector sensitivity by applying a peak flux cut at 0.4 ph sec -1 
cm~ 2 leaving a total sample of 555 LGRBs. 

Figure U (lower panel) shows the best fitting results to 
the Swift differential peak flux distribution. The peak flux 
intervals, number of bursts and burst rates data used for 
the fit is given in Table IA2I We find best fit results of po — 
0.09iHi G P c "V"\£» = S.Ol^oa x 10 52 ergs _1 anda = 
-3.8±o;g with a goodness of fit y^jdof = 19.9/10 ~ 2. 
T hese values are in agreement w ith the no-evolution model 
of ISalvaterra fc Chincarinil (|2007l ). 



Cmax is the rate in the second brightest detector. We finally 
select long bursts with T90 > 2s and with set a threshold of 
fmin = 0.4 corresponding to the value at which the trigger - 
ing algorithm is almost 100% efficient (|Paciesas et al.|[l999tl . 
We find that 596 GRBs meet these criteria. 

Table 17.11 summarizes the different data segments and 
the total number of bursts meeting the above constraints. 
We show the trigger numbers (rather than the dates) and the 
total detector run time to form each segment. In our later 
analysis we will apply the log P - log T method to both the 
combined data and to the largest of the individual segments. 

Figure|4](top panel) shows the best fitting results to the 
BATSE differential peak flux distribution using Q = 0.677T 
and a live time of AT = 3.19 yr. The peak flux inter- 
vals, number of bursts and burst rates data used for the 
fit are given in Table IA1I We find a best fit of po — 



0.12tnniGpc- 3 yr- 1 ,L, = 3.1±f* x 10 51 ergs _1 and a = 



^taoi^pc "yr = H.lLiy x iu ei 

-2 2^ 017 ^t-*" ~rfij- ^ ^,2 1 



-32 



. The goodness of fit is given as \ jdof = 6.6/13 



0.51. 



7.2 



The Swift log N — log P distribution 

In recent studies Izhand (|201ll ) and iBromberg et all (|2012l ) 
have suggested that the T90 = 2 s division of short and 
long GRBs based on t he BATSE bimodial distribution 
|Kouveliotou et ai1 [T993) is a detector dependent categori- 
sation and therefore not appropriate for Swift bursts. Other 
studies have suggested an interme diate duration class o f 
bursts between these two classes (|Horvath et al.l l2010bl ). 
For our LGRB sample, we obtain peak flux values from 
the Swift online catalogue^ but rather than employing a 
T90 cut, we use the catagorisations given in the Jochen 
Greiner online catalogue (JG) of well localized GRBs0. As 
the burst catagorisations in this catalogue are subject to 
review through follow up observations we find it a use- 
ful resource to obtain our LGRB sample. We find that 
from the Swift peak flux sample of 649 bursts up to burst 
120224A, 644 are catagorised in the JG catalogue. We ex- 
clude the three SL-GRBs, 060218, 060505 and 100316D 
(see section 110.11 for further discussion of these bursts) 
and 12 bur sts catagorised as SGRB-EE (a ll but one have 
T 90 > 2s) (|Norris et al.ll201ll ). We follow iGuetta fc Piranl 



http: / /swif t .gsf c.nasa. gov/docs/swif t/archive/grb_table/ 



THE LOG P 
LGRBS 



LOG T DISTRIBUTION OF 



In this section we apply the estimated parameters, L*, a 
and po to the 90% PEH log P - log T thresholds given in 
section UTTl and attempt to constrain peak flux PEH data. 

As a distant bright burst could produce a similar peak 
flux as a closer burst of moderate brightness, an intrinsic 
scatter will be present in P(T) data. This, in combination 
with the spatial distribution of bursts, means it is difficult 
to use a 90% log P- log T threshold to separate burst pop- 
ulations. The spatial dependence introduced in section 16.21 
will provide a stronger case for untangling burst popula- 
tions. The goal of this section is two fold: a) to show that 
the parameters obtained using the log N - log P fits are com- 
patible and complementary with the log P- log T relation; 
b) to demonstrate how data conditioning methods improve 
the quality and quantity of the PEH data. 
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Figure 5. The top panel shows the BATSE P(T) sample obtained 
by chronologically combining the data sets A-F of table P7~T1 The 
bottom panel shows the same data after applying the FromMAX 
algorithm to temporally re-ordering the data from the maximum 
peak flux. The PEH data are indicated in each panel by squares. 
The top panel illustrates how an energetic event occurring at 
early observation-time can decreases the PEH sample. The bot- 
tom panel shows how some simple processing using FromMAX 
can circumvent this problem. 
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12 3 

Log T (days) 

Figure 6. The logP-logT data of Figure[5]is constrained using 
cosmological 90% PEH curves. At high observation-time the data 
is also constrained by a Euclidean model. 



8.1 The BATSE log P — log T distribution 

The log P- log T method can be used to fit to either the 
individual data runs of BATSE shown in Table 17.11 or a 
combined set - this allows several consistency tests of the 
method. We firstly use all the available data to construct a 
single time series P(T) by chronologically combining the 6 
data sets A-F. A temporal dependence is obtained by adding 
the observation-times of the jth data set Tj to the maximum 
of the previous max(Tj_i). When joining two consecutive 
data sets an additional factor of T = 1.5 days, the mean 
time between events for BATSE, is also added to the first 
event of Tj to ensure a sufficient delay time is included. 

Figure [5] shows the combined BATSE log P(T) data. 
The top panel shows the pre-processed sample with PEH 
data shown by squares. We see that as observation-time in- 
creases, so too does the probability of a large event. The plot 
also illustrates how a bright event can occur near the start 
of a P(T) time-series (at around 10 days) and thus minimise 
the PEH sample - we will see later that this effect can also 
be seen in the individual data sets A-F. 

The bottom panel of Figure [5] shows the same data after 
applying the fromMAX technique to temporally re-order the 
data from the brightest event (as described in section^. The 
advantage of the fromMAX method is clearly illustrated - 
the PEH sample is increased from 5 events to 14 events. 

To test that PEH data set obtained through From- 
MAX (top panel) is statistically compatible with that ob- 
tained from the pre-conditioned data (top panel), we ap- 
ply a 2D-Kolmogorov-Smirn ov test (2DKS) (|Peacocklll983l ; 
iFasano fc Franceschinilll987f ) and obtain a KS probability of 
29% showing the two samples are statistically consistent. 

Figure [6] shows the 90% log P- log T thresholds and 
peak flux PEH data of BATSE. We see that the data is well 
constrained by the cosmological threshold. In addition to the 
increased size of the sample, the improvement obtained us- 
ing the FromMAX method can be further validated through 
binomial maximum likelyhood (BML) estimates for obtain- 
ing data within the 90% thresholds. We find a BML of 92% 
for the data shown in Figure [6] (1 failure in 14). We find a 
BML of 60% for the pre-processed PEH data shown in the 
top panel of Figure [5] (2 failures in 5). We note that a bias 
towards the lower threshold (e = 95%) is evident - this is a 
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Figure 7. As for Figure [5] but this time showing on the left 
PEH data from pre-processed P(T) samples B,C,F and G. The 
FromMAX PEH data from each sample is shown on the right. We 
see that the sample size is increased in all but the first case. 



result of the data conditioning rather than an intrinsic bias, 
as will become apparent when we next apply the method to 
individual data sets. 

Figure [7] shows a repeated analysis for the 4 largest 
individual data segments, B,C,F and G. The left sided panels 
display the pre-processed data (PEH data is again shown by 
squares) and the right sided panels show the time series after 
applying the FromMAX method. We see that the PEH data 
has increased in 3 of the 4 samples. 

Figure [8] shows the 90% log P - log T thresholds applied 
to the 4 sets of PEH data shown in the right hand panels of 
Figure [7| The plot shows that the PEH data is again well 
constrained. The bias towards the lower threshold (e = 95%) 
evident in Fig. is only apparent for data set F, suggesting 
that this is a manifestation of the data conditioning rather 
than the BATSE sample. The plot does suggest that meth- 
ods su ch as data sp l itting and recombining procedures as 
used in I Howell et al.l l|2007al ) on simulated data could prove 
useful to increase PEH data sets from long time-series. 

Table [2] compares the amount of PEH data and the 
BML estimations obtained using the 4 sets of pre-processed 
data with those obtained after applying the FromMAX 
method. All cases show that the amount of PEH data 
and/or the BML estimates have improved from the process- 
ing. Additionally, the 2DKS probabilities are given, which 
all indicate good statistical compatibilities between the pre- 
processed PEH samples and after applying FromMAX. 



8.2 The Swift log P — log T distribution 

Figure [9] shows the Swift FromMAX time-series - the pre- 
processed P(T) data was shown previously in Fig.[T] In com- 
parison to BATSE sample shown in the top panel of Fig. fSJ 
one could argue that the Swift data already seems well con- 
ditioned to extract PEH data. This is because: a) the largest 
event in the the pre-processed Swift data occurred near the 
end of the observation-time; b) no significantly bright event 
has occurred at early observation-time. 

Figure [TO] shows the 90% log P - log T thresholds us- 
ing both the pre-processed and FromMAX time-series. Both 
PEH samples recorded BML values of 100% although the 
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Table 2. A comparison of PEH data obtained from pre- 
conditioned data with that obtained using the FromMAX 
method. The table shows that the amount of data and/or bi- 
nomial maximum likelyhood (BML) estimates imrove with the 
FromMAX method. A 2D Kolmogorov-Smirnov test on the two 
PEH samples from each data set shows that the statistical PEH 
signature is not lost by using the FromMAX method. 



Data set 




B 


C 


F 


G 


Pre-conditioned 
PEH data 


Size 
BML 


5 

80% 


5 

100% 


4 

100% 


4 

75% 


FromMAX 
PEH data 


Size 
BML 


5 

100% 


8 

100% 


8 
88% 


6 

100% 


2DKS Probability 




0.84 


0.85 


0.68 


0.28 




12 3 

Log T (days) 

Figure 8. As for Figure \E\ but with the PEH data of Figure 
[7] obtained using FromMAX. We see that the PEH data is well 
constrained by the 90% log P - log T thresholds. 



FromMAX method has produced an increase in data. The 
statistical compatibility between the PEH samples before 
and after applying FromMAX is confirmed by a value of 
0.46. 

We have shown in this section that a log P- log T dis- 
tribution of GRBs is in agreement with a log N — log P dis- 
tribution and that the method can be applied to data sets 
obtained from different detectors. We have further shown 




10 10 10 

observation time (days) 

Figure 9. The P(T) distribution of Swift after applying the 
FromMAX method. 



how some processing can improve the quality and quantity 
of data. In the next section we shall apply the same princi- 
ples to the redshift domain and will show how the technique 
can be used to differentiate between different populations of 
burst. 



9 THE SWIFT LOG Z- LOG T DISTRIBUTION 

As discussed earlier, the PEH sample of a Z(T) distribution 
appro aches the 1ow-jz regime rapidly (|Coward &: Burmanl 
2005), thus the log Z-log T method is most sensitive to 
the closest occurring events. The dependence on only the 
closest events mean that selectio n effects such as the 'red- 
shift dessert' |Coward et alJlioOct ) can be ignored. Addition- 
ally, as the source rate evolution is reasonably well predicted 
within z ~ 2 the choice of SFR model will not change the 
curves significantly. For the closest events, as redshift de- 
terminations would be expected to be more accurate, we 
include values obtained through all methods: absorbtion, 
emission and photometrically. 

Figure [TT] illustrates how the high-z variations have 
minimal influence on the log Z- log T curves by examin- 
ing the redshift selection function. The top panel shows 
the function N(< Zl) defined by equation [16] assuming 
Pl = 0.4 - one may recall that log Z — log T curves are pro- 
duced by introducing a temporal dependence (equation [TT]). 
The dashed line shows the same function without the factor 
S(L) = /^ ma3 |p L &(L)dL which can be referred to as a de 



tecto r dependent scaling factor (see eq. 1 of ICoward et al.l 
2008). The plot shows that this function has little effect until 
z ~ 3 (z ~ 2 for the BATSE instrument) - as the detection 
threshold decreases this value of z increases. 

The lower panel shows the log Z-log T for these two sce- 
narios showing that there is little change in the two curves. 
Therefore one can make reasonable estimates without an 
accurate description of the luminosity distribution of the 
sources. For a log Z-log T function, the rate density is the 
dominant variable. As such, one can discriminate between 
transient populations of different rate density solely by their 
spacial dependence. 
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Figure 10. The 90% log P- log T thresholds for Swift are used 
to constrain PEH data. Top panel - PEH data from the pre- 
processed time-series. Bottom panel - PEH data using the From- 
MAX time-series. 
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Figure 11. The top panel shows the function N(< Zl) with 
and without the scaling factor discussed in the text. The lower 
panel shows that the two log Z — log T curves differ little and are 
dependent only on the closest events of a source population. 
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Figure 12. The Swift Z(T) distribution with PEH data indicated 
by squares. The top panel shows the pre-processed sample. The 
lower panel shows the PEH data has increased after applying the 
FromMIN procedure. 



Figure shows the Swift Z(T) distribution with PEH 
data indicated by squares. The Swift Z(T) sample of 173 
bursts is obtained by taking bursts from the selection of 
section \T7Z\ with certain redshift measurements. The data is 
given in Table IB1I The top panel shows the pre-processed 
data and the lower panel shows the data after applying 
the FromMIN procedure described in section [4] The Figure 
shows that the PEH data set has been increased by applying 
FromMIN (the PEH data is indicated by bold in Table iBTj) . 

Figure [T3l shows that the PEH data is well constrained 
by the 90% log Z-log T thresholds. To account for the Swift 
efficiency in obtaining a redshift, po is scaled by a factor 
of 0.3 to account for the number of bursts with redshifts 
(172) from our overall sample (576). To test the compatibil- 
ity of the PEH data obtained using the FromMIN method 
with that obtained from the pre-processed Z(T) we perform 
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Figure 13. The 90% log Z - log T thresholds and PEH data 
for the Swift LGRB sample obtained using the FromMIN 
method. The data is well constrained using the rate po = 



0.09 



+0.01 
-0.01 



Gpc 3 yr 1 estimated through a log N - log P fit. 



a 2DKS test. A probability of 0.87 indicates that the two 
samples are from the same distribution. 

We have shown in this section by extracting PEH data 
from the Swift Z(T) sample that the log Z- log T relation is 
in good agreement with both the log N— log P and logP- 
logT relations. We have also illustrated that a log Z — log T 
distribution relies strongly on the spacial distribution of the 
sources and is invariant to the detector selection function. 
In the next section we will demonstrate how the method can 
differentiate between source populations. 



10 THE LOG Z - LOG T RELATION APPLIED 
TO SL-GRBS 

10.1 The log Z — log T relation as a prob of burst 
populations 

In this section we will demonstrate that the log Z - log T 
method can be used to differentiate between different source 
populations. A key principle we will use is that PEH data 
represents the rarer, brighter or closer events that are less 
subject to detector selection bias. Thus, insight into the av- 
erage properties of a source population can be provided from 
only a small sample of PEH events. 

Bursts from the same distribution, with the same in- 
trinsic rate density, should produce PEH data that is con- 
strained by a log Z- log T threshold. A separate burst pop- 
ulation, occurring at a higher rate, will produce rarer PEH 
events at earlier observation-times and visa-versa for pop- 
ulations at lower rates. For two mixed Z(T) populations, 
their observation-time sequence will differ, thus producing 
outlying events to a 90% log Z — log T threshold. 



10.2 The Swift sub-luminous GRB sample 

In addition to the two excepted SL-GRBs in the Swift 
sample, GRB 060218 and GRB 100316DA, a number of 
other candidates have also been discussed. Based on their 
low lumi nosities and subsequent Po isson detection prob- 
abilities, Wan derman fe Piranl l)2010l l suggested three ad- 
ditional bursts could contribute to the SLGRB popula- 
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GRB 


^peak 


Tgo 


T 




a 


z 




phs _1 cm 2 


s 


days 


ergs -1 






060218 


0.25 


2100 


427 


1.4 x 10 47 


-2.3 


0.03 


060505 


2.65 


4 


504 


9.1 x 10 48 


-1.3 


0.089 


100316D 


0.1 


1300 


1915 


4.8 x 10 47 


-2.3 


0.059 


Table 3. 


Data for the 


three 


Swift sub-luminous 


bursts. 


The 



isotropic rest frame (1-10000 keV) luminosities are obtained from 
the peak luminosity values (ph s _1 ) using a simple power law fit 
with values of a taken from the Swift online catalogue. 



tion: GRB 050724, GRB 051109 (see also iBromberg et all 
l2012f l and GRB 060505. Of these three bursts we include 
GRB 060505 in our sample. Although GRB 050724 was a 
long duration burst (Tgo ~ 96 s), further a nalysis has sug- 
gested that it is of the SGRB category JCampana et aD 
120061 ; iGrupe et al.ll200d; iMalesani et al.ll2007l )- most likely 
a SGRB-EE |Norris et alj|20ld . l201ll ) GRB 051109 is also 
omitted due to uncertainty in the redshift estimatior0 (in 
accordance with our selection criteria of section I12|) . 

Despite GRB 060505 being relatively nearby and well 
observed, intense photometric and spectroscopic searches 
found no evidence of an associated supernova despite the 
fact an event ~ 100 times fainter than SN 1998bw wou ld have 
been detected (|McBreen et al.l 120081 : feu et all 120091'). Dust 
obscu ratio n was excluded fo r this burst and iFvnbo et al.l 
(2006) and lOfek et all (|2007l ) argued that the simplest ex- 
planation was that GRB 060505 was simply the closest 
observed SGRB. Other studies however found the proper- 
ties o f the host galaxy consistent with that of the LGRB 
class l|Ofek et al.ll2007l:lThoene fc Fynbdl2007l ; iThone et aJj 
120081 ). Furthermore. iMcBreen et alj (|2008l ) showed the spec- 
tral lag was consistent with that of a 

lgreFL 

and an 

extens ive study of the afterglow performed by Ixu et al.l 
(2009) found afterglow parameters within the range for other 
LGRBs. In view of its luminosity, we include this burst from 
our analysis of SL-GRBs in the next section. However, to 
further address the classification of this burst, we will inves- 
tigate its compatibility to the Swift SGRB population in the 
final section. Table [3] shows the main parameters for these 
three SL-GRBs which will be used in the analysis of this 
section. 



10.3 The log Z — log T analysis of Swift SL-GRBs 

Figure IT41 shows the P(T) time-series for the LGRB sample 
obtained using FromMIN and including the three SL-GRBs 
060218, 060505 and 100316D. The PEH data are indicated 
by squares - the SL-GRB sample are clearly shown as sig- 
nificantly close events. On inspection, GRB 060505 stands 
out as a particularly rare event due to its occurrence at a 
relatively early observation-time. 

Figure [15] shows the 90% log Z - log T thresholds based 
on the LGRB parameters used in section [9] along with the 
PEH data. This plot clearly shows that the three SL-GRBs 



0.08 was reported in GCN 5387 



9 A putative host galaxy at z 

by D. A. Perley et. al. 

^ LGRBs support the lagluminosity relation of iNorris et all 

feOOd ) while SGRBs have zero lag. 
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Figure 14. The Z(T) time series for Swift LGRBs including three 
Swift SL-GRBs. 
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Figure 15. The log Z-log T distribution for Swift LGRBs includ- 
ing three SL-GRBs which are labeled. The dashed line represents 
a 0.02% probability for obtaining a LGRB within redshift 2 at an 
observation-time T. 



are outliers to the distribution suggesting they are from a 
population occurring at a higher rate. Based on the 90% 
thresholds of Fig. 1131 assuming that the SL-GRBs are from 
the low-probability tail of the LGRB distribution, one would 
expect to wait 87 (1800) years to observe a burst as close as 
GRB 060505 (GRB 060218). The probability of observing 
these bursts at time T can be estimated by extending the 
PEH threshold upper limit through e - this new limit is 
shown in Fig. [15] by the dashed curve. We find a probability 
value of e = 0.00015 is required to constrain the SL-GRB 
bursts at the LGRB rate. T his value is consistent with that 
estimated bv lCowardl |2005l ) for the probability of observing 
the low-z burst GRB 980425. 

If the SL-GRB sample are from a distinct population 
one should be able to constrain them using seperate 90% 
log Z- log T thresholds set at a higher rate value po.SL and 
a LF representative of the SL-GRB sample. For the LF 
we consider the form used for the LGRB sample but with 
L* = 10 49 erg s _1 which is around three orders of mag- 
nitude lower than that derived in section [7] As described 
in the previous section, consideration of S(L) is not essen- 
tial but is included for completeness. To estimate po.SL for 
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Figure 16. The log Z-log T distribution for Swift SL-GRBs. 



GRB 


Peak Flux 


redshift 


101219A 


4.1 


0.718 


100206A 


1.4 


0.4068 


100117A 


2.9 


0.92 


090510 


26.3 


0.903 


070724A 


2.03 


0.457 


061217 


2.36 


0.827 


051221A 


12 


0.5465 


050509B 


3.71 


0.226 


080905At 


6.03 


0.1218 


060505* 


2.65 


0.089 



Table 4. The Swift SGRB data sample. The lower two entries 
are investigated as members of the SGRB population: * possible 
SL-GRB but no SN observed to deep limits; < low-z outlier to the 
Yonetoku relation with high energy properties typical of a SGRB 
(see section 10.4). 



the S w jft SL-GRB sample we follow iGuetta fc Delia Valid 
lj2007h : fCoward et alJ l|2012h and calculate 



P0,SL = ^ 



l l l 



(19) 



Here V ma x is the maximum volume each burst could be de- 
tected, T is the maximum observation-time for the sample, 
Q is the sky coverage and we assume r\ z =0.3 as for the 
LGRB sample. To determine 7 max we use equation 0] and 
assume a flux limit of 0.15 ph sec - - this corresponds with 
the detection threshold of 95% of the Swift LGRB sam- 
ple. We obtain p ,sL = 147ta9°Gpc~ 3 yr~ 1 where the er - 
rors are the 90% Poisson confidence limits (Gehrels 1986). 
This rate, based solely on the Swift sample, is in agree- 
ment with ot her studies of the SL- GRB population (e.g. 
Coward! 120051; ISoderberg et alJlioOrj; IGuetta fc Delia Valid 
20071; iPian et al.ll2006l ; iLiang et al.ll2007l : I Virgilii et alj|200a 
Howell et all 120 lift 

Figure [16] shows the 90% log Z - log T thresholds cor- 
responding to po.SL- To show that the log Z — log T method 
allows constraints to be made solely on po.SL, the dashed 
lines show the curves obtained without the selection function 
S(L). The SL-GRB PEH data is well constrained demon- 
strating good agreement with the rate estimates and sup- 
porting the hypothesis of a distinct SL-GRB population of 
bursts. Additionally, the figure supports the hypothesis that 
GRB 060505 is a member of the SL-LGRB category. 



10.4 Testing the connection between SL-GRBs 
and XRFs 

We note that the two sub-luminous bursts GRB 060218 
and GRB 100316D are often categorised as XRFs. Al- 
though XRFs are generally understood to r epresent the 
fainter and softer part of the GRB distribution (|Zhang||2007l ; 
ISakamoto et al.l l2008h the relation between XRFs and SL- 
GRBs is not clear. Looking at Swift data the durations of 
both GRB 060218 (2100 s) and GRB 060218 (1300 s) are 
both much higher than those of other XRFs in the cat- 



alogut[^] suggesting that a different mechanism separates 
the two categories. It is therefore interesting to extend the 
log Z — log T analysis of the previous section to investigate 
the connection between SL-GRBs and XRFs using Z(T) 
data. We find that including two XRFs with secure red- 
shiftiQ XRF 050416A (z=0.65) and XRF 050824 (z=0.83); 
does not change the result given in Figure 1151 Our anal- 
ysis using the Swift sample with secure redshifts therefore 
suggests that XRFs and SL-GRBs are from a different pop- 
ulation. 



10.5 A log Z — log T analysis of GRB 060505 

Although the last section shows that based on a log Z - log T 
analysis, GRB 060505 is not of the LGRB category, the lack 
of an accompanying SN to stringent limits still suggests it 
could be of the SGRB class. A additional test is therefore to 
include this event in the Z(T) distribution of Swift SGRBs 
and repeat the previous analysis. 

To calculate the observed rate of SGRBs, po.s, we add 
GRB 060505 to the SGRB sample of lCoward etall (|2012T ) 
and extend equation 1191 through the factor Rb/s = 6.7 to 
account for the reduced sensitivity of Swift for detecting 
SGRBs in comparison to BATSE. We omit an anomalous 
burst GRB 080905A from our sample - a repeat analysis 
at the end of this section will consider this burst. Setting 
rj z = 9/41 we obtain p ,s = 8.571 7 4 7 7 Gpc- 3 yr -1 - 

Figure[T7]shows the 90% log Z - log T thresholds includ- 
ing GRB 060505, which is shown as a clear outlier. We find 
that a value of e = 0.007 is required to constrain this burst 
indicating a > 99% probability that a SGRB would not have 
occurred at z = 0.09. Our analysis therefore suggests that 
GRB 060505 belongs in the SL-GRB sample of bursts. 

The previous analysis excluded the burst GRB 
080905A. This spectrally hard, short duration burst (T90 = 
1 s) showed negligible spectral lag, no evidence of extended 



11 Of the 11 XRFs in the Swift data up to January 2012 the 
longest T90 is 103s; 6 have durations below 10s. 

12 As in section [7. 21 we again base our catatgorisations and red- 
shifts on the Jochen Greiner online catalogue. 
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Figure 17. The log Z- log T distribution for Swift SGRBs show- 
ing that GRB 060505 is an outlier to this distribution. The dashed 
curve indicates a > 99% probability that a SGRB would not have 
occurred at z = 0.09. 




2 3 4 5 

Log T (days) 

Figure 18. As for Fig.fTflbut including the burst GRB 080905A. 
We see that this burst, as well as GRB 060505 are outliers. The 
dashed curve indicates a ~ 1% probability that GRB 080905A 
could have occurred at z = 0.12. 



emission or of an associated SN to deep limits - all prop- 
erties of a SGRB resulting from a compact object merger. 
It was spectroscopically associated with a galaxy at red- 
shift z = 0.1218 making it the close st possible SGRB 
l|Rowlinson et al.ll2010l ; iKann et alJl201ll L 

Figure [18] shows the result of a repeated analysis of 
GRB 060505 including GRB 080905A with 90% \ogZ- 
logT thresholds corresponding to a recalculated rate of 
Po,s = 13.5^12 2Gpc -3 yr _1 . The plot shows that as well 
as GRB 060505, GRB 080905A is also an outlier to 
the 90% thresholds. This finding supports a n ana lysis by 
iGruber fc for the Fermi /GBM collaboration! (|2012h which 
found GRB 080905 A to be a clear outlier to the Yone- 
tok u, Ep.rest-Lp relation llYonetoku et al.l |2004 ); but see 



also iBorgonovo fc Biornssonl ( 20061 ). They found that GRB 
080905A would require a redshift z ~ 0.9 to be consistent 
with the Yonetoku relation, suggesting that the host, a large 
almost face-on s p iral g alaxy, could be a foreground galaxy. 
iRowlinson et al.l (|201dh noted that the offset of the after- 
glow was large (18.5 kpc) but comparable to other SGRB 
locations, especially considering the relative size of the host. 



Our analysis suggests a ~ 1% probability (indicated 
by the dashed curve) that GRB 080905A could have oc- 
curred at z = 0.12. Furthermore, in agreeme nt with 
IGruber fc for the Fermi/GBM collaborationl (|2012l ). we find 
that given a redshift of z ~ 0.9, this burst would be con- 
strained by the log Z - log T relation. 

Given that the observed prompt properties of this burst 
are highly consistent with a SGRB and the redshift is rea- 
sonable secure, the above analysis is interesting and implies 
that this burst may deserve additional attention. It certainly 
further highlights the ability of the log Z- log T method to 
select out bursts with anomalous properties from a given 
sample. 



11 CONCLUSIONS 

In this study we have shown how the redshift and peak-flux 
distributions of GRBs have an observation-time dependence 
and have derived two new relations from the standard inte- 
gral distributions - the log P - log T and the log Z - log T. We 
have shown how a PEH data set can be extracted from the 
redshift and peak-flux time series and constrained by 90% 
thresholds using the above two relations. We have demon- 
strated that the FromMIN and FromMAX methods improve 
the PEH data sample and circumvent a bias introduced by 
the occurrence early rare event. 

By applying the log P - log T method to both Swift and 
BATSE data we showed that the relation is in agreement 
with parameters obtained through a log N — log P fit. We 
then showed that the same parameters could be used to 
constrain PEH data using a log Z- log T relation. By in- 
cluding the Swift SL-GRB sample of bursts we showed that 
this sample were outliers to the LGRB sample and could 
be constrained by a 90% log Z — log T threshold at a higher 
rate. We suggest that this is further evidence that SL-GRBs 
are a discrete population of bursts. 

We have shown that the methods presented here are 
strongly dependent on the event rate density of the sources 
and a good description of the scaling function S(L) is not 
essential. Additionally, in the redshift regime where the 
method is dependent on only the closest events, high-z biases 
do not effect the analysis. Used in the peak flux regime, the 
method depends on only the brightest of bursts which are 
less likely to be missed because of detection thresholds. We 
suggest that an observation-time analysis is a useful com- 
plement to other methods that use the whole distribution 
and can be a good indicator of selection bias. 

The log P- log T method, in comparison with a logZ - 
logT technique can not be used to discriminate between 
populations occurring at different rates. A close burst could 
arise from the low luminosity tail of the LF and equivalently 
a distant burst would be intrinsically bright. Therefore a 
log Z — log T relation which has only a spacial dependence is 
a stronger indicator. As an intrinsic luminosity is determined 
from both the redshift and peak flux, a log L — log T relation, 
which includes a measure of both the energy and spatial dis- 
tribution could prove useful. Such a relation may be useful to 
untangle populations such as SGRBs and SGRB-EEs which 
have similar redshift distributions. 

We analysed the burst GRB 060505 using the logZ- 
logT relation to determine if it could be of the SL-GRB 
category. Including this burst with the other two Swift SL- 
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GRBs, GRB 060218 and GRB 100316D, we found that all 
three bursts were constrained at rate po.SL = 147Gpc _3 yr _1 
estimated independently by a Vmax analysis. Therefore, in 
addition to its low-luminosity and relative close proximity a 
log Z - log T analysis supports the hypothesis that this burst 
is a SL-GRB. In contrast to other SL-GRBs, despite intense 
photometric and spectroscopic searches, this burst showed 
no evidence of an associated supernova. To tests its compat- 
ibility with the SGRB population, we performed a further 
log Z — log T analysis and found a greater than 99% proba- 
bility that a SGRB would not have occurred as close as this 
burst. 

If the progenitor of GRB 060505 was not a com- 
pact object merger, upper l imits on the e j ected 56 Ni 
mass of M( 56 Ni)~ 1O~ 3 M (|Tominaga et al.l 120071 ) sup- 
port the possibility that this burst may have been accom- 
panied by a low energy SN undergoing sig nificant fallback 
l|Wooslevfc Weaverll 19951 ; iFrver et aljfcood ) . 

An additional analysis of the burst GRB 080905A 
shows how the log Z- log T relation can select out bursts 
with anomalous properties from a given sample. Using the 
method in a bootstrapping type scheme may be a useful 
consistency test for both the completeness of sample se- 
lections and some of the spectral energ y relations related 
to the observed quantities of GRBs (e.g. Amati et al.l |2002| ; 
lYonetoku et al. 1 12004 iGhirlanda et al.ll2005f) . 
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APPENDIX A: THE DIFFERENTIAL PEAK 
FLUX DISTRIBUTION DATA OF BATSE AND 
SWIFT LGRBS 



PI 

phsec -1 cm~ 2 


P2 

ph sec -1 cm 


N 

- 2 yr ~i 


AN/AP 
yr _1 ph sec cm 2 


0.4 


0.51 


25 


222 


0.51 


0.66 


22 


145 


0.66 


0.84 


20 


113 


0.84 


1.1 


17 


88.2 


1.1 


1.4 


15 


61 


1.4 


1.8 


13 


38.4 


1.8 


2.3 


11 


24.6 


2.3 


2.9 


8.8 


15 


2.9 


3.7 


7.2 


8.78 


3.7 


6.1 


11 


4.64 


6.1 


7.8 


3.8 


2.44 


7.8 


10 


3 


1.5 


10 


13 


2.1 


1.06 


13 


16 


1.9 


0.415 


16 


21 


1.5 


0.324 


21 


11 


2.7 


0.129 


Table Al. The 


data used in 


fitting the differential peak fh 


distribution of the BATSE long GRB sample 


PI 

ph sec -1 cm~ 2 


P2 N 
phsec - 1 cm -2 yr — 1 


AN/AP 
yr — 1 phsec cm 2 


0.4 


0.56 


9.7 


59 


0.56 


0.8 


7.8 


33.4 


0.8 


1.1 


11 


34.6 


1.1 


1.6 


8.9 


19.1 


1.6 


2.2 


11 


16.3 


2.2 


3.2 


6.7 


7.2 


3.2 


4.5 


3.8 


2.87 


4.5 


6.3 


3.6 


1.95 


6.3 


8.9 


.3.8 


1.44 


8.9 


13 


1.8 


0.49 


13 


18 


0.97 


0.187 


18 


25 


1.7 


0.227 


25 


50 


0.97 


0.0388 



Table A2. The data used in fitting the differential peak flux 
distribution of the Swift long GRB sample 



APPENDIX B: PEH REDSHIFT DATA 
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Table Bl. The Swift redshift-observation-time data Z(T) used in the analysis. The data is taken from 
December 19 2004 and re-ordered as a time-series using the FromMIN method. The PEH data is shown as 
bold. 



GRB 


^obs 


rcdshift 


GRB 


T u 

1 obs 


rcdshift 


GRB 


- 1 obs 


rcdshift 


050904 


8 


6.29 


071117 


811 


1.331 


091018 


1512 


0.971 


050908 


12 


3.344 


071122 


816 


1.14 


091020 


1514 


1.71 


050922C 


26 


2.198 


080129 


888 


4.349 


091024 


1518 


1.092 


051016B 


50 


0.9364 


080207 


896 


2.2 


091029 


1523 


2.752 


051109A 


73 


2.346 


080210 


899 


2.641 


091109A 


1533 


3.076 


051111 


75 


1.55 


080310 


929 


2.42 


091127 


1551 


0.49 


060115 


144 


3.53 


080319B 


938 


0.937 


091208B 


1562 


1.063 


060124 


153 


2.296 


080319C 


938 


1.95 


100219A 


1638 


4.6667 


060202 


161 


0.783 


080330 


949 


1.51 


100302A 


1651 


4.813 


060206 


165 


4.048 


080411 


960 


1.03 


100316B 


1665 


1.18 


060210 


169 


3.91 


080413A 


962 


2.433 


100418A 


1697 


0.6235 


060223A 


182 


4.41 


080413B 


962 


1.1 


100425A 


1704 


1.755 


060418 


237 


1.489 


080430 


979 


0.767 


100513A 


1722 


4.772 


060502A 


251 


1.51 


080520 


999 


1.545 


100621A 


1760 


0.542 


060510B 


259 


4.9 


080603B 


1012 


2.69 


100724A 


1793 


1.288 


060512 


261 


2.1 


080604 


1013 


1.416 


100728B 


1797 


2.106 


060522 


271 


5.11 


080605 


1014 


1.6398 


100814A 


1813 


1.44 


060526 


275 


3.221 


080607 


1016 


3.036 


100901A 


1830 


1.408 


060604 


283 


2.68 


080707 


1046 


1.23 


100906A 


1835 


1.727 


060605 


284 


3.78 


080710 


1049 


0.845 


101219B 


1938 


0.55 


060607A 


286 


3.082 


080721 


1060 


2.591 


110106B 


1960 


0.618 


060707 


316 


3.425 


080804 


1073 


2.2045 


110128A 


1982 


2.339 


060708 


317 


1.92 


080805 


1074 


1.505 


110205A 


1989 


2.22 


060714 


323 


2.711 


080810 


1079 


3.35 


110213A 


1997 


1.46 


060729 


338 


0.54 


080905B 


1104 


2.374 


110213B 


1997 


1.083 


060814 


353 


0.84 


080906 


1105 


2.1 


110422A 


2066 


1.77 


060904B 


373 


0.703 


080913 


1112 


6.695 


110503A 


2077 


1.613 


060906 


375 


3.686 


080916A 


1115 


0.689 


110715A 


2149 


0.82 


060908 


377 


1.8836 


080928 


1127 


1.692 


110731A 


2165 


2.83 


060912A 


381 


0.937 


081007 


1136 


0.5295 


110801A 


2165 


1.858 


060926 


395 


3.2 


081008 


1137 


1.9685 


110808A 


2172 


1.348 


060927 


396 


5.47 


081028A 


1157 


3.038 


110818A 


2182 


3.36 


061007 


406 


1.261 


081029 


1158 


3.8479 


111008A 


2232 


4.9898 


061021 


420 


0.3463 


081118 


1177 


2.58 


111107A 


2261 


2.893 


061110A 


439 


0.758 


081121 


1180 


2.512 


111209A 


2293 


0.677 


061110B 


439 


3.44 


081203A 


1192 


2.05 


111228A 


2312 


0.714 


061121 


450 


1.314 


081222 


1211 


2.77 


111229A 


2313 


1.3805 


061126 


455 


1.1588 


081228 


1217 


3.4 


120119A 


2338 


1.728 


061222A 


481 


2.088 


081230 


1219 


2 


050126 


2381 


1.29 


061222B 


481 


3.355 


090102 


1226 


1.547 


050223 


2408 


0.5915 


070110 


504 


2.352 


090205 


1259 


4.6497 


050315 


2430 


1.949 


070208 


532 


1.165 


090313 


1297 


3.375 


050318 


2433 


1.44 


070306 


560 


1.4959 


090417B 


1331 


0.345 


050319 


2434 


3.24 


070318 


572 


0.836 


090418A 


1332 


1.608 


050401 


2446 


2.9 


070411 


595 


2.954 


090423 


1337 


8.26 


050505 


2480 


4.27 


070419A 


603 


0.97 


090424 


1338 


0.544 


050525A 


2500 


0.606 


070506 


620 


2.31 


U9U426 


1340 


2.609 


050603 


2508 


2.821 


070521 


635 


1.35 


090429B 


1343 


9.2 


050730 


2565 


3.967 


070529 


643 


2.4996 


090516A 


1360 


4.109 


050801 


2566 


1.56 


070611 


655 


2.04 


090519 


1363 


3.85 


050814 


2579 


5.3 


070612A 


656 


0.617 


090529 


1373 


2.625 


050820A 


2585 


2.612 


070721B 


695 


3.626 


090530 


1374 


1.3 


050826 


2591 


0.297 


070802 


706 


2.45 


090618 


1392 


0.54 








070810A 


714 


2.17 


090715B 


1419 


3 








071003 


767 


1.6043 


090726 


1430 


2.71 








071010A 


774 


0.98 


090809 


1443 


2.737 








071010B 


774 


0.947 


090812 


1446 


2.452 








071020 


784 


2.145 


090814A 


1448 


0.696 








071031 


795 


2.692 


090926B 


1490 


1.24 








071112C 


806 


0.823 


090927 


1491 


1.37 
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